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Abstract 

We study well-posedness and asymptotic dynamics of a coupled sys- 
tem consisting of linearized 3D Navier-Stokes equations in a bounded 
domain and a classical (nonlinear) full von Karman shallow shell equa- 
tions that accounts for both transversal and lateral displacements on a 
flexible part of the boundary. We also take into account rotational in- 
ertia of filaments of the shell. Out main result shows that the problem 
generates a semiflow in an appropriate phase space. The regularity 
provided by viscous dissipation in the fluid allows us to consider simul- 
taneously both cases of presence inertia in the lateral displacements 
and its absence. Our second result states the existence of a compact 
global attractor for this semiflow in the case of presence of (rotational) 
damping in the transversal component and a particular structure of 
external forces. 

Keywords: Fluid-structure interaction, linearized 3D Navier-Stokes 
equations, nonlinear shell, global attractor. 
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1 Introduction 



We consider a coupled (hybrid) system which describes an interaction of 
a homogeneous viscous incompressible fluid which occupies a domain O 
bounded by the (solid) walls of the container S and a horizontal (flat) bound- 
ary f2 on which a thin (nonlinear) elastic shell is placed. The motion of the 
fluid is described by linearized 3D Navier-Stokes equations. To describe 
deformations of the shell we use the full von Karman shallow shell model 
which accounts for both transversal and in-plane displacements. For details 
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concerning the shell model chosen we refer to [421 |21>1 EZ] and also to the 
papers [281 [291 (SB [331 EHl [36] and the references therein. 

This fluid-structure interaction model assumes that large deflections of 
the shell produce small effect on the fluid. This corresponds to the case 
when the fluid fills the container which is large in comparison with the size 
of the plate. 

We note that the mathematical studies of the problem of fluid-structure 
interaction in the case of viscous fluids and elastic plates/bodies have a long 
history. We refer to [SI [I3l [191 [20l EH [24] and the references therein for 
the case of plates/membranes, to [15] in the case of moving elastic bodies, 
and to [H [U O El [?1 [16] in the case of elastic bodies with the fixed interface; 
see also the literature cited in these papers. 

We also note that the global (asymptotic) dynamics in nonlinear plate- 
fluid models were studied before in [91 H3] . The article [9] deals with a class 
of fluid-plate interaction problems, when the plate, occupying f2, oscillates in 
longitudinal directions only. This kind of models arises in the study of blood 
flows in large arteries (see, e.g., [19] and the references therein). A fluid- 
plate interaction model, accounting for only transversal displacement of the 
plate, was studied in [13J. In contrast our mathematical model formulated 
below takes into account both transversal and in-plane displacements. 

Let O C R 3 be a bounded domain with a sufficiently smooth boundary 
dO. We assume that dO = $7 U S, where 

fl C {x = (zi; x 2 ; 0) : x' = (x 1 ;x 2 ) G R 2 } 

with a smooth contour T = d£l and S is a surface which lies in the subspace 
_ | X3 < q| 'j 1 ^ ex t er i or normal on dO is denoted by n. We have that 
n = (0;0; 1) on O. 

We consider the following linear Navier-Stokes equations in O for the 
fluid velocity field v = v(x,t) = (v 1 (x,t);v 2 (x,t);v 3 (x,t)) and for the pres- 
sure p(x, t): 

v t -uAv + Vp = G f (t) in Ox(0,+oo), (1) 

divu = in Ox(0,+oo), (2) 

where v > is the dynamical viscosity and Gf(t) is a volume force (which 
may depend on t). We supplement (P) and ([2]) with the (non-slip) boundary 
conditions imposed on the velocity field v = v(x,t): 

v = on S; v = (v ; v 2 ; t> 3 ) = (u t ; u 2 ; wt) on £1, (3) 

where u = u(x,t) = (u 1 ;u 2 ;w)(x,t) is the displacement of the shell occu- 
pying $7. Here w stands for transversal displacement, u = (u l ;u 2 ) — for 
lateral (in-plane) displacements. 

To describe the shell motion we use the full von Karman model which 
takes into account the rotational inertia of the filaments and possible pres- 
ence of in-plane acceleration terms (see the literature cited above). 
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We denote by Tf(v) the surface force exerted by the fluid on the shell, 
which is equal to Tn\n, where n is a outer unit normal to dO at O and 
T = {Tij}f j =1 is the stress tensor of the fluid, 

T i:j = Tij (v) = v (v l x . + vQ - p5ij , i, j = 1, 2, 3. 
Since n = (0; 0; 1) on Q, we have that 

T f (v) = (K<+<);Ki+<);2<f 3 -p). 

Below for some simplification we assume that for the case considered Young's 
modulus E and Poisson's ratio u G (0, 1/2) are such that Eh = 2(1 + fi), 
where h is the thickness of the shell. In this case the (elastic) stress tensor 
{Nij} is given by the formulas 

2 2 
Nu = YZTJl ( £n + , N 22 = Y^Jl ( £22 + ^ £ ^) ' ^12 = £ 12, (4) 

where the deformation tensor {£ij} has the form 

£ n = «ii + k i w + ^(^i) 2 ' 

£22 = u 2 X2 + A; 2 w + ^(^a) 2 , 
£12 = ul 2 +ul 1 +w Xl w X2 . 

Here k\ and A;2 are curvatures of the initial form of the shell which are 
sufficiently smooth functions of x' £ Q. 

After an appropriate rescaling of the parameters and functions we can 
model shell dynamics by the following equations 

M a (w tt + jw t ) + A 2 w + hN u + k 2 N 22 

- d Xl (Nuw Xl + N 12 w X2 ) - d X2 (N 12 w Xl + N 22 w X2 ) 

= G 3 (t)-2i/^t) 3 +p in Ox(0,oo), (5) 

where M a = 1 — aA, and 

Qui = 9 X1 N U + 8 X2 N 12 + Gi(i) - u(vl 3 +v xi ), 

gui = d Xl N 12 + 8 X2 N 22 + G 2 (t) - v{v 2 xz + v X2 ), (6) 

where G s h{t) = (G\]G 2 ;Gz){t) is a given body force applied to the shell. 
a > and g > are constants which take into account rotational inertia 
and in-plane inertia of the shell, respectively, 7 is a non-negative parameter 
which describes intensity of the viscous damping of the shell material. 
We impose the clamped boundary conditions on the shell 



1 1 2 1 & w 

u \dQ = u \ d Q = W\ d Q = — 

on 



= (7) 

an 
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and supply (H])-© with initial data for the velocity field v = (v 1 ; v 2 ; v 3 ) and 
the shell displacement vector u = (u 1 ; u 2 ; w) of the forrd3 

vo, u\ t=Q = u , w t \ t=0 = w 1 , g [ut\ t=0 -ui] = 0, (8) 



t=0 



where u = (u ,u ). Here vq = (vq;v 2 ;vq), uq = {uq\u^;wq), w\, and u\ 
(u\;u 2 ) are given vector functions which we specify later. 

We note that ([2]) and ([3]) imply the following compatibility condition 
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w t (x',t)dx' = for all t > 0. (9) 
This condition fulfills when 

/ w(x',t)dx' = const for all t > 
Jn 

and can be interpreted as preservation of the volume of the fluid. 

We emphasize that even in the linear case we cannot split system ([!])- 
(jU) into two sets of equations describing longitudinal and transversal plate 
movements separately, i.e., we cannot reduce the model considered to the 
cases studied in [13]. The point is that the surface force Tf(v) is not 
the sum of the corresponding loads in the models [9] and [13j . The models 
in [HI [13] are much simpler in several respects. For instance, in the case of 
longitudinal plate deformations only (see [9]) the equations which correspond 
to ([6]) do not contain the terms v x . and the model does not require any 
compatibility conditions like ([9]) because the volume of the fluid obviously 
preserves in the case of longitudinal deformations. In the case of purely 
transversal displacements |13] the force exerted on the plate by the fluid 
contains the pressure only. 

In the following remark we describe some structural properties of the 
shell model chosen. 

Remark 1.1 (A) One can see that the equations for the in-plane displace- 
ment vector u = (tr; u 2 ) can be written in the vector form as follows: 

Qu tt + Au = B(w) + (Gi(t) - v(vl a + v Xl );G 2 (t) - u(v 2 X3 + v 3 X2 )), (10) 

where the operator A in f)10[) is defined by 

/(l + \)d 2 1+ d 2 X2 Xd xlX2 \ 

^d XlX2 d 2 Xl + (1 + \)d 2 J 



A 



with A = i±g and D(A) = [H 2 (Q) f| H^(n)] 2 . The nonlinear term B{w) in 
(1101) has the form 



^ l _ d Xl \2xxw + (w Xl ) 2 + /j,(w X2 ) 2 ] + d X2 [w Xl w X2 ]^ 
B(w)= ^ J 

d Xl [w Xl w X2 ] + _ d X2 [2x 2 w + (w X2 ) 2 + fj,(w Xl ) 2 ] 



1 We put the multiplier g in the fourth relation of © to emphasize that this relation 
is not needed in the case of negligibly small in-plane inertia (g — 0). 
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where x\ = fci+yuA^ and X2 = kz+fiki. Thus for fixed w and v the equations 
for the in-plain displacement u are similar to the standard equations of 2D 
(linear) elasticity theory. 

(B) The equations in ([5]) and © can be also written in the form 

(1 - aA)(w tt + 7^t) + A 2 w + trace {K,jV{u)} 

= div {jV{u)Vw} + G 3 (t) - 2vd X3 v 3 + p 

and 

Qu tt = div 



where K = diag (k±, hi) and 

/ N u N 12 \ 

with u = (ui; 112), ^(e) = 2(1 — [/i trace e • / + (1 — /i)e] and 
eo(^) = ^(Vn + V T n), /(s) = i S ®s, s£l 2 . 

This form of the full von Karman system was used earlier by many authors 
in the case when the fluid velocity field v is absent and k{ = (see, e. g., [33] 
or [23] and the references therein). 

Our main goal in this paper is to prove well-posedness of the problem 
in (HI)-© in the class of finite energy solutions (see Theorem I3.3|) and to 
show a possibility of compact long-time dynamics (see Theorem 15.11 on the 
existence of a compact global attractor). We note for the proof of unique- 
ness in Theorem 13.31 relies substantially on the i? 1 -regularity of wt, which 
follows from the structure of the mass operator M a and involves Sedenko's 
method (see |35} [36J ) . To prove the existence of a global attractor in The- 
orem ED we use J. Ball's method (see [5] and [32]). To apply this method 
we need the property 7 > 0, i.e., assume a presence of rotational damping 
in the transversal component of displacement. The question whether the 
system under consideration demonstrates compact long-time behavior with- 
out mechanical damping in the shell component is still open. In contrast 
we note that the existence of global attractors in the models considered in 
[13] does not require any mechanical damping and compact asymptotic 
dynamics of the corresponding system is guaranteed by viscous dissipation 
of fluid. One of the reasons for this is that the models in [91 [13] do not 
include higher order (rotational type) inertial terms and thus the finiteness 
of the full dissipation integral for the displacement follows from the viscosity 
of the fluid via the compatibility condition in ([3]) . 

The paper organized as follows. In Section[2]we provide some preliminary 
material related to Sobolev spaces and the Stokes problem. In Section [3] we 
state and prove our main well-posedness result. Section [5] deals with long- 
time dynamics. 
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2 Preliminaries 



In this section we introduce Sobolev type spaces we need and provide with 
some results concerning the Stokes problem. 

2.1 Spaces and notations 

To introduce Sobolev spaces we follow approach presented in [40] . 

Let D be a sufficiently smooth domain and s£R. We denote by H S (D) 
the Sobolev space of order s on the set D which we define as a restriction 
(in the sense of distributions) of the space H s (W l ) (introduced via Fourier 
transform). We define the norm in H S (D) by the relation 

\\u\\l jD = inf jlMlsRd : weH s (R d ), w = u on L>j 

We also use the notation || • ||o = || • ||o,d and (-,-)_d for the corresponding 
L2 norm and inner product. We denote by Hq(D) the closure of C^°(D) in 
H S (D) (with respect to |[ • \\ s ,d) and introduce the spaces 

Ht{D) := [f\ D : / G H s (R d ), supp/ C d} , s G R. 

Below we need them to describe boundary traces on 17 C dO. We endow 
the classes H*(D) with the induced norms ||/||* D = \\f\\ s R d for / G H^(D). 
It is clear that 

\\f\\ s ,D < \\f\\* s ,n, f&m(D). 

However, in general the norms || • \\ St D and || • ||* D are not equivalent. It is 
known that (see [ffil Theorem 4.3.2/1]) that C^\D) is dense in H*(D) and 

Ht(D) C H^(D) C H S {D), s G M; 
H S (D) = H S (D), -oo<s<l/2; 

Ht{D) = H°(D), -l/2<s<oo, s- 1/2 ^{0,1,2,...}. 

In particular, H*(D) = H$(D) = H S {D) for \s\ < 1/2. Note that in the 
notations of [31] the space H™ +1 ^ 2 (D) is the same as H^ +1 ^ 2 (D) for every 
m = 0, 1, 2, . . . , and for s = m + a with < a < 1 we have 




where d{x,dD) is the distance between x and dD. The norm || • ||* D is 
equivalent to || • \\ s ,d in the case when s > — 1/2 and s — 1/2 {0, 1, 2, . . .}. 

Understanding adjoint spaces with respect to duality between C^°(D) 
and [Cq°(D)]' by Theorems 4.8.1 and 4.8.2 from 00] we also have that 

[m(D)}' = H- S (D), s G M, and [H 8 (D)]' = H~ S (D), s G (-oo, 1/2). 
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Below we also use the factor-spaces H S (D)/M. with the naturally induced 
norm. 

To describe fluid velocity fields we introduce the following spaces. 

Let ^{O) be the class of C°° vector- valued solenoidal (i.e., divergence- 
free) functions on O which vanish in a neighborhood of S. We denote by X 
the closure of ^(O) with respect to the L2-norm and by V the closure with 
respect to the i? 1 (0)-norm. One can see that 

X = {v = (v 1 ;^ 2 ;^ 3 ) G [L 2 {0)f : div v = 0, j n v = (v, n) = on S} ; 

(11) 

and 

V = {v = {v l ;v 2 -v z ) G [H l (0)f : div v = 0, v = on S} . 

We equip X with L2-type norm || • \\q and denote by (•, -)o the corresponding 
inner product. The space V is endowed with the norm || • ||y = ||V • For 
some details concerning this type spaces we refer to [39], for instance. 

We also need the Sobolev spaces consisting of functions with zero average 
on the domain O, namely we consider the space 

L a (n) = g L 2 (n) : J u {x')dx' = o| 

and also H S (Q) = H s (£l) n ^2(^) for s > with the standard F s (fi)-norm. 
The notations Hl{Vt) and Hq{£1) have a similar meaning. 
To describe shell displacement we use the spaces 

W = Hl{VL) x Hl(VL) x flg(fi), y = L 2 (0) x L 2 (fi) x (12) 

where i? a equals to -ffo(^) w ith the equivalent norm ||-||# = |Hln+ a ll^*lln 
and corresponding inner product. 

Remark 2.1 Below we also use Hq (O) as a state space for the displacement 
of the plate. It is clear that Hq(VL) is a closed subspace of Hq(CI). We denote 
by P the projection on Hq(Q) in iJg(O) which is orthogonal with respect 
to the inner product (A-,A-)q. One can see that (/ — P)Hq(Q) consists of 
functions u G Hq(Q) such that A 2 u = const and thus has dimension one. 

2.2 Stokes problem 

In further considerations we need some regularity properties of the terms 
responsible for fluid-plate interaction. To this end we consider the following 
Stokes problem 

—vAv + Vp = g, div v = in O; 
v = on 5; v = tp = (ip 1 ; ip 2 ; V> 3 ) on fl, (13) 
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where g G [L 2 (0)] 3 and ip G [L 2 (Q)] 2 x L 2 (f2) are given. This type of bound- 
ary value problems for the Stokes equation was studied by many authors 
(see, e.g., [25] and [39] and references therein). We collect some properties 
of solutions to f|13|) in the following assertion. 



Proposition 2.2 The following statements hold. 

(1) Let g G [H- 1+ °(0)f, and i() G [hI /2+0 (ft)] 3 with J n ip 3 (x')dx' = 0. 
Then for every < a < 1 problem fjl3[) /tos a unique solution {v;p} 
in [H l+u {0)f x [H a {0)/R] such that 

\\ v \\[Hi+"(0)} 3 + \\p\\h°(0)/R < CO {llffll[J3'-i+^(0)]3 + IHI[^+i/2 (c) ] 3 } ■ 

(14) 

(2) lfg = 0,ifj£ [H~ l,2+a (tt)f , < <r < 1, / n V 3 (x')^x' = 0, tten 

IMI[H CT (0)P + lbllH-i+-(C)/R < co||^|| [i? -i/s +CT(n)]3 . (15) 

In particular, we can define a linear operator Nq : [L 2 {Q)\ 2 x L 2 (f2) ^ 
[i? 1 / 2 ^)] 3 6y £/te formula 

Ar , f -z/Aio + Vp = 0, divu; = in O; , , 

N ip = w iff < _ „ . „ (16) 

( to = onS; !» = i|) on!l, v ' 

for ip G [L 2 (Q)} 2 x L 2 (ft) (2Vo^ solves JT3J) toift g = Oj. Jt /oZZows 
/rom Jj^[ ) and $15]) that 

N : [H s ^n)} 2 x H°(n) h-> [^(Ojfnl continuously 

for every —1/2 < s < 3/2. 

Proof. We just combine proof of Proposition 2.1 [9] and Proposition 2.2 
|13j . The argument in these references rely on the consideration in [251 [39] 
and also in [18]. □ 

3 Well-Posedness Theorem 

To define weak (variational) solutions to ([I])-® we need the following class 
Ct of test functions on O: 



<P G L 2 (0,T; [H^O)] 3 ), <p t G L 2 (0,T; [L 2 (C)] 3 ), 
div0 = O, <f)\s = 0, ( p\ n =b=(b 1 ;b 2 ;d), 
d G L 2 (0,T;£ 2 (ft)), 6* G L 2 (0, T; j = 1,2, 

dt G L 2 (0,T;H^(n)), ti> t G L 2 (0, T; L 2 (0)), j = 1,2. 



We also denote £^ = {<f> G £ T : 0(T) = 0}. 
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Definition 3.1 A pair of vector functions (v(t);u(t)) with v = (v 1 ;?; 2 ;?; 3 ) 
and u = (u 1 ; u 2 ; u>) is said to be a weak solution to the problem in {I])-© 
on a time interval [0, T] if 

. veL oo (0,T;X)r\L 2 (0,T;V); 

• u G L^T-Hliti) x x tf 2 (O)); 

• ut £ L2(0,T; \h1^ 2 (Q)] 3 ) and the compatibility condition = 
(«|; u 2 ; w t )(t) holds for almost all i G [0, T]; 

• (gu\;gu 2 ;w t ) G L oo (0,T;L 2 (^) x L 2 (fi) x H a ) and u(0) = u ; 



for every G £y with 0[q = b = (b ; b 2 ; d) the following equality holds: 



r 



-/ (v,<j) t )odt + v E(v,4>)odt- I [{M a w t ,dt)o, + g(u t ,b t )n]dt 
Jo Jo Jo 

+ 7/ (M a w t ,d) n dt + (Aw,Ad) n dt+ a{u,b) n dt+i q{u,b) n dt 
Jo Jo Jo Jo 

= (v , <j)(0)) o + {M aWl ,d{0)) n + g(ui, 6(0)) n 

+ [ T (G f (t),<f ) ) dt+ [ T (G sh (t),b) Q dt, (17) 



o 

where u = (u ; u 2 ), 6 = (ir; 6 2 ) and also 
1 3 

=2 E (<+<■) (4+4,), 
2 l + 



a(«,6) = VVW, V6% + — ^(divu, div6) n , (18) 

' 1 — u 

1=1 r 

g(«,6) =(^1^11+^2^22,^)^ 

+ {Nww Xl + iV^^jdaJn + (N 12 w Xl + N 2 2W X2 ,d X2 )n 

+ t— ; [(KJ 2 + ^K) 2 . &ijn + (K 2 ) 2 + mKi) 2 , On] 
i /i 

2 

+ (m^m^b^ +& 2 Jn + — — (tO.Klft^ +Ki6 2 2 ) n - 

1 /i 

Remark 3.2 (1) In the case when we neglect the inertia of longitudinal 
deformations (g = 0) the equations in (|6|) (or in (jTDjl ) become elliptic. How- 
ever we keep the initial data for the in-plane displacement (u l \u 2 ). The 
point is that the first order evolution for (u 1 ;?/ 2 ) goes from the boundary 
condition for the fluid velocity in ([3]). 

(2) It also follows from ([3]) and from the standard trace theorem that 
for every weak solution (v(t);u(t)) we have that 

\Ht)\\l l/2{n) + iKWII^ + lk 2 (*)llij/. (n) < C\\Vv(t)\\ 2 (19) 
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for almost all t G [0, T]. This estimate does not depend on g > and provide 
us an additional regularity estimate for in-plane shell velocities even in the 
case g = 0. Below we use this observation to suggest unified way to prove 
well-posedness result not distinguishing cases g > and g = (in contrast 
with [42] , see also [35] (g = 0) and [36] (e > 0)). 

(3) Taking in (fTT)) </>(£) = J* t T x{ T )dr ■ tj), where x is a smooth scalar 
function and ?/> belongs to the space 

V = {V G F | = /3 = (/3 1 ; /3 2 ; 5) € i^) x x H$(n) } , (20) 

one see that the weak solution (v(t);u(t)) satisfies the relation 

{v(t),ijj)o + (M a w t (t),5)n + g(u t (t)Jh 



[vE{v, V) + (Aw, A<5) n + a(u, p) + j(M a w t , S) n 

+ q(u,p)-(G f ,i/>)o-(G ah ,p)n]dT (21) 

for all t G [0,T] and ip G V" with = £ = (P^P^S) and /3 = (/3 1 ;/? 2 ). 

As a phase space we use 

( {(vq;uo;ui) e X x W x Y : vo = ui on £1} , g > 0, 
% = < r „ 1 (22) 

I {(«o;no;wi) G X x W x H a : (u ) = wi on f!} £ = 

with the norm 



||(vo;uo;ui)||^ = ||« ||o + \\uo\\w + IkilliL + dK^L^i" 1 " ' 



f2- 



where and y are given by (IT21) . We also denote by % a subspace in % of 
the form 

ft = {{vo-uo-m) en-. w e # 2 (^)} , 

where wq is the third component of the initial displacement vector uq. 
Our main result in this section is the following well-posedness theorem. 

Theorem 3.3 Assume that U = (vq\Uq;ui) G H, G f (t) G L 2 (0,T;V), 
G sh (t) G L 2 (0,T; [H- l / 2 {9)] 2 x H' 1 ^)), 7 > 0. We also assume that 
g > (in the case g = the data U\ are not fixed). Then for any interval 
[0,T] there exists a unique weak solution (v(t);u(t)) to {I])-® with the 
initial data Uq. This solution possesses the following properties: 

• In the case g > we have 

U(t;U Q ) = U{t) = {v{t)-u(t)-u t (t)) G C(0,T;X xWxY), (23) 
If g = 0, then 

U{t;U ) = U(t) = (v(ty,u(t);w t (t)) G C(0,T;X xWx H a ). (24) 
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The solutions depends continuously (both in strong and weak sense) on 
initial data. More precisely, if g > and U n — > Uq in the norm of % 
(resp. weakly in %), then U(t;U n ) — > U(t;Uo) strongly (resp. weakly) 
in % for each t > 0. In the case g = the corresponding convergence 
take place in X x W x H a . 

The energy balance equality 



£{v{t),u{t),u t {t)) + v E(v,v)dT + j \\w t \\ 2 Ha dT 
Jo Jo 

= £(v ,u , til) + / (Gf,v)odT+ / (G sh ,u t )ndT (25) 
Jo Jo 

is valid for every t > 0, where the energy functional £ is defined by the 
relation 



£(v,u,u t ) 



1 



\l + Wl^wtWl + g\\u t \\l + ||AH|n + Q{i 



(26) 



with 
Q{u) 



7 1^M) In ( £ " + + 2/X£ll£22 + \ <yl ~ 

= 2(TT^) L ^ + N<22 ~ 2 ^ NllN22 + 2(1 + ■ (27) 

Proof of Theorem EOl 

We start with the following elliptic type property of the functional Q(u). 

Proposition 3.4 There exists a positive constant C such that for every 
u = (u l ]u 2 ]w) & W = Hl(9) x Hl(9) x H$(Q,) we have that 



Q(u) + \\w\\ 2 n + \\w u4 



3/2, n 



(28) 



Proof. One can see that 



[<] 2 <24 + 2 



kiw + ^ (w x J 2 



1 2 



, [ul 2 + u 2 Xl ] 2 < 2e 2 12 + 2 [w Xl w X2 Y 



Therefore using the embedding H 1 / 2 ^!) C L 4 (0) we obtain that 



/ f«) 2 + «) 2 + K 2 + <; 

Jn L 



dx' < C 



Q{u) + \\w\\n + IMl3/2,n 



Thus, property (|28p follows from the Korn inequality (see, e.g., Theorem 3.1 
in Chapter 3 of [IT]). □ 
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Now we use the compactness method and split the argument into several 
steps. 

Step 1. Existence of an approximate solution for the case g > 0. For the 
construction of Galerkin's approximations in this case we use the same idea 
as in [13J (which was inspired by the method developed in [8] for the case of 
a linear plate interacting with nonlinear Navier-Stokes equations). 

Let {V'ilieN be the orthonormal basis in X = {v G X : i;|„ = 0} consist- 
ing of the eigenvectors of the Stokes problem: 

-At/ji + Vpi = ^iipi in O, divipi = 0, ipi\do = 0. 

Here < fix < fj,2 < ■ • ■ are the corresponding eigenvalues. Denote by 
felieN the basis in Hq(CI) which consists of the eigenfunctions of the fol- 
lowing problem 

(A&, Aw) n = ki(M a ^,w) n , Vw G tf 2 (ft), 

with the eigenvalues < k\ < ki < . . . and such that (M a £j,£j)Q = Sij. 
Further, let {?7i}ieN be the basis in Hq(TI) x Hq(Q) which consists of eigen- 
functions of the problem 

a(r]i,w) = ki(rn,w)a, Vw G Hq(£1) x Hq(U), 

with the eigenvalues < k\ < ki < • • • and ||?7i||n = 1> The form 0,(17, w) ^ s 
given by (fT8j) . 

Let fa = -/Vo(0; 0; and 0j = Nofa; 0), where the operator ./V"o is defined 
by (fl6|) . One can see that ??j G [^T 3 / 2 ~ 5 (fi)] 2 . Therefore by Proposition I2.2I 
we have that fa, fa G [F 2 ~ 5 (C)] 3 n F for every 5 > 0. 

In what follows we suppose fa = fa, Q = (77^ ; 0) , Ki = ki, and 7$ = 0; 
4> n+i = fa, Cn+i = (0,0, K n+i = Kf, and 7„ +i = 7 for i = 1, . . . , n. We 
also define the parameter £>, by the relations: Qi = Q for & = 1, . . . ,n and 
= 1 for i = n + 1, . . . , In. 

We define an approximate solution as a pair of functions (v nm ; u n ): 

m 2n 

(t) = ^a i (i)^ + ^/3 7 (t)<^, 
i=l j=i 

2n 

= X)/9i(<)Ci + (0;0; (/ - P)w ), (29) 
which satisfy the relations 

2n 2n 

+5^/5,-(t)(0j,^ fe )o + v/j, k a k (t) + $j(t)E(<j)j,ip k )o = (Gf,Tp k )o 
i=i i=i 

(30) 
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for k = l...m, and 

m 



2n 



i=l 



2» 



+ v ati(t)E(il>i, Mo + »Y1 Pj i<t>k)o + 

i=l j=l 

+ lkPk{t) + q(u n {t), Ck) = (G f (t), 4> k )o + (G sh (t), ( k )n (31) 

for k = 1, . . . , 2n. This system of ordinary differential equations is endowed 
with the initial data 

(0) = U m (v - N (0; 0; Wl )) + N (0; 0; P nWl ), 

u n (Q) = (R n (ul; Uq); P n Pw + (I - P)w ), u n (0) = (R n (u\;uj); P n w x ), 

where H m is the orthoprojector on Lin{vpj : j = 1, ...,m, } in X, P n is 
orthoprojector on Lin{£i : i = 1, . . . , n} in L2{Q) and i? n is orthoprojector 
on Lin{r]i : i = l,...,n} in £2^) x L2(fJ). Since II m , i? n and P n are 
spectral projectors we have that 

(«v,m(0);«n(0);wn(0)) («o;«o;«i) strongly in "ft as m, n oo. (32) 
We can rewrite system ([30]) and (f3Tj) as 



dt\m) 



+ g(a(t),/3(t)J(t)) = G(t) 



for some locally Lipschitz function g : flj m + 4n i_ >• ]R" 1+2ra , where the function 
G lies in L 2 (0, T; R m+2n ) and 



M 



"0 0" 




St 


+ 



{(^,^)o}r fc =i mAk)o)7£=i 
m^ih}tk=i {(Mi))o}f k=1 



(33) 



where M = diagj^i, . . . , £?2n}- The first matrix in ([33]) is nonnegative and the 
second one is symmetric and strictly positive (since the functions : 
i = l,...,m,j = l,...,2n} are linearly independent). Therefore system 
(f30j) and (j3Tj) has a unique solution on some time interval [0, T'\. 
It follows from ([29]) that 



'Jn,rn 



i=i 



where Nq is given by ()16j) . This implies the following boundary compatibility 
condition 

v n ,m(t) = d t u n {t) on O. (34) 
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Step 2. A priori estimate and limit transition, g > 0. Multiplying (|30|) 
by cfcfc(t) and f)31 j) by f3k(t), after summation we obtain an energy relation of 
the form ([25]) for the approximate solutions (v ntm ;u n ) (for a similar argu- 
ment we refer to [13] and also to the classical source [12] on (non-interacting) 
shell evolution). By Proposition 13.41 and relation (|19h this implies the fol- 
lowing a priori estimate: 

SUp lll^m (t)\\h+ 
te[o,T] 



sup 

te[o,T] 



\M l J 2 d t W n {t)\\l + Q\\d t U n {t)\\l + ||Au> n (t)||n + ll«n 



+ ^ ||V«„, m (t)|&dt + jf \\d t u n (t)\\ 2 [Hl/2m2 dt<C T (35) 

for any existence interval [0, T] of approximate solutions, where the constant 
Ct does not depend on n and m. In particular, this implies that any ap- 
proximate solution can be extended on any time interval by the standard 
procedure, i.e., the solution is global. It also follows from (|35|) that the 
sequence {(v n ^ m ; u n ; dtu n )} contains a subsequence such that 

;u n ;d t u n ) — (v; u; dfU) * -weakly in L^iO, T; 7i), (36) 
weakly in L 2 (0,T;V). (37) 

Moreover, by the Aubin-Dubinsky theorem (see, e.g., [37, Corollary 4]) we 
can assert that 

u n -»• u strongly in C(0, T; H*~ e (n) x H^- £ (n) x H 2 - e (n)) (38) 
for every e > 0. Besides, the standard trace theorem yields 



H^ 2 (dO)^ , (39) 



u n ,m - 1 v weakly in L 2 10, T; 
thus, we have 

d t u n <9 t u weakly in L 2 (0, T; Hl ,2 {ti)). (40) 
One can see that (fn,m! 

;dtu n )(t) satisfies (fl~7l) with the test function 

(j) of the form 

p q q 

<t> = <h,q = Z}7i(*M + X) +J2vj(t)$j, (41) 

i=l j=l j=l 

where p < m, q < n and ji, 6j, rjj are scalar absolutely continuous functions 
on [0, T] with time derivatives from £2(0, T), such that ji(T) = 5j(T) = 
rjj(T) = 0. Using ([36]) and ([37|) . we can easily pass to the limit in linear 
terms. Limit transition in the nonlinear terms can be done the same way as 
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in [12] (see also [33] or [23J for the same type argument). Thus, one can show 
that (v;u;dtu)(t) satisfies (fT7|) with (ft = 4> Ptq , where p and q are arbitrary. 
By (j32j) and (|38|) we have u(0) = uq. The compatibility condition ([3]) follows 
from flM}, dUJI and p9l) . 

We conclude the proof of the existence of weak solutions by showing that 
any function (f> in Ct can be approximated by a sequence of functions of the 
form (I4ip . we refer to [13] for a similar argument in the case of zero in-plane 
deformations. 

This solution satisfies the energy inequality: 

£(v(t),u(t),u t (t)) + v I \\Vv\\ldT + 1 l \\w t (r)\f Ha dT 
Jo Jo 

t r t 



< £{vq,uq,ui) + / (G f (r),v)odT + / {G sh {r), u t )udT 



for almost all t > 0, where the energy functional 8 is defined by (|26p . 

Step 3. Existence of weak solutions in the case g = 0. We fix some 
u% from L,2(p.) and consider the corresponding solution (v e (t);u g (t)) with 
g > 0. As above, Proposition 13.41 and relation (I19p imply the following a 
priori estimate: 



ess sup 

te[o,T] L 



MOIIo + l|My2^(t)||2 + \\A We (t)\\n + \\u e (t)\\U 

+ £ \\Vv e (t)f Q dt + £ \\dtn e m 2 [Hl /2 m2 dt < C T (42) 

for any interval [0,T], where the constant Ct does not depend g € (0,1). 
This implies that the family {(v g ; u e ; dtu g )} contains a subsequence such 
that 

(v g )Ug) — (v)u) * -weakly in L oo (0,T;X x W), 
d t w Q -± d t w * -weakly in Loo(0, T; Hq(Q)), 
v g — v weakly in 1^(0, T; V), 

when g — > 0, and also 

u e ->> u weakly in L 2 (0, T; # 1/2 (SO) ) , 

a t n e d t u weakly in L 2 (o,T; H^ 2 (dO) ^ 

Therefore we have that 

u e -> u strongly in C(0,T; H^~ £ {Q) x i? 1_£ (O) x #o~ e (^)) 

for every e > 0. This allows us to make limit transition when g — > in (|17p 
and prove the existence of weak solutions for the case g = 0. 
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Step 4- Uniqueness. We use the same idea as in [3B]. However the 
additional regularity of ut (see Remark 13.2( 2)) makes it possible to cover 
both cases g > and g = simultaneously. 

Let (v(t);u(t)) and (v(t);u(t)) be two solutions to the problem in ques- 
tion with the same initial data. These solutions satisfy (J42]) and their dif- 
ference (v(t);u(t)) = (v(t) — v(t);u(t) — u(t)) possesses properties 

v = (v 1 ; v 2 ; v 3 ) G .MO, T; X) f] L 2 (0, T; V); 

u = [u x \u 2 \w) G (0, T; H 1 (Q) x x flg(J2)); 

^ G Loo(0, T; Hq,), K 1 ;^ 2 ) G L 2 (0, T; [^ 1/2 (0)] 2 ); 

and by (|2ip satisfies the relation 



(v(s),ip)o + (M a w t (s),5)n + g(ut(s),P)n 

= - [ [vE(v, rp) + (Ate, A<5)n + a(S, ft) 
Jo 

+ 7(M a w t , <y) n + ?(«, /3) - g(tt, /?)] (43) 
for all s G [0,T] and V G V with = Z 3 = (/3 1 ;/3 2 ;5) and /3 = OS 1 ;/? 2 ), 



where V is defined by ([20]) , 

One can see that u(t) lies in Hq(JV) x i?o(f2) x -^o(^) f° r eacn * — an< ^ 

■0 = i[i>](s) = / v(a)da G V" with tp\ = ft = u(s) for each s G [0,T], 
•/() " ; 

where we have introduced the notation 

f s 



i[h](s) = ( h(a)da, s G [0,T], 
Jo 



for any integrable function h(t) with values in some Banach space. 

Substituting this t/j in (I43p for the fixed s G [0, T] after integration from 
to t over s we obtain that 

£(*) + 0l|fi(t)lln<2 t ds [ S dT[q(u(T),u(s))-q(u(T),u(s))} . (44) 
JO JO 

where 

C(i)=IKH(i)ll^ + l|Mi/Vi)|| 2 , + 2 i / / '(fa£?(*[t,]( fl ),iH(«))) 

j 

+ 2 7 / ||MV 2 ^ s )|| 2 d S + ||AiH(i)||^ + «(i[n](t),i[u](i)). 

JO 

Now we estimate the integral term in (I44p . We first note that the term 
q(u(r), u(s)) has the structure 

4 

g(«(r),«(s)) = J^^(«(r), «(«)). 
i=i 

Here 
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• q±(u, u(s)) is a linear combination of terms of the form (x tB, w(s))q 
and (xnDiu i ,w(s))ci, where D/ = d Xl and xq, xu are smooth function, 
i,* = 1,2; 

• ^(w, u(s)) is a linear combination of terms of the form 
(D l wD i w,w(s)) n , (DiwDiwD k w,D m w(s)) n , (xi m wD t w, D m w(s))n; 

• qs(u,u(s)) is a combination of the terms (D^w, Diw, D m u l (s))n with 
2 = 1,2; 

• (74(2,14(5)) consists of the terms (D^u 1 , Diw, D m w(s))n with i = 1,2. 
We also denote 

I m {t) = [ ds f d,T[q m (u(T),u(s)) -q m (u(T),u(s))] 
Jo Jo 

d,T[q m (u(T),i[u](t) -i[u](r)) - q m {u(T), i[u)(t) -i[u](r)]. 





Below the notation a ~ 6j means that a is a linear combination of terms 6j 
and a < bi means that a can be estimated by a linear combination of hi. 
With these notations we have 

h(t) ~ / ds / [(x w(r),u;(s))c) + {x H Diu j {t) , w(s))n] 
Jo Jo 

= / (is [(x i[u>](s), w(s))n) + (mDii[u j ](s),w(s))n\ . 
Jo 

This implies that 

\h(t)\<C f\(s)ds, te[0,T]. (45) 
Jo 

The most complicated term in (72 is 

r(to,«;(s)) := (DiwDiwD k w, D m w(s))n- 

We consider 

I[r] = j ds j dr [r(w(T),w(s)) — r(w(r),w(s))] 
Jo Jo 

~ / dT(Aw*(r)Aw*(r)I>fc«;(r),I> m i[«;](t)-D m i[«;](T)) n , 

JO 

where is either w) or w. Since i/ 1 (Q) C L P (Q) for every 1 < p < 00 

|/[r]|<C /*[||w(r)||l > n + ||fi?(r)||l >n ] 
Jo 

x |k(r)|| 1)Q [|KH(i)|| 2)Q + |KH(r)|| 2i n]dT 

<e||iM(t)ll2,n + C £ f [\\w(r)\\l n + ||zH(r)||| n ] dr 
Jo 
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fore every e > 0, where the constant C e depends on T and on bounds for 
solutions w and w. 

Using a similar argument in other terms of 6/2 we obtain the estimate 

\h(t)\ < e€(t) + C e [* €(s)ds, t€[0,T\, Ve>0. (46) 
J o 

To estimate I$(t) we note that 

h(t) ~ /V (Ato*(r)2?fc«;(r),Z) m *[^](t) - Z) m i[^'](r)) n , 



o 

where w* is either w or u; as above. To estimate 13(f), we use the following 
Bresis-Gallouet type inequality 

\\fg\\n < Cl {log(l + a)} 112 ||/||Q|bl|i,n + r^-H/lli.nllfflli.n (47) 

1 + 0" 

for every f,g€ H 1 ^) and for all cr > (this inequality can be proved the 
same way as Lemma A. 3. 6 |12j . see also the Appendix in [14]). 
Inequality (|47p and a priori bound (|42j) imply 

||L> z w*AcHh <ci {log(l + cr)} 1/2 ||^*||2,nlkl|i,n + 7-^— ||u>*||2,n||Hl2,fi 

1/2 1L..11 , C 2 (T) 



Thus, 



< Ci(T) {log(l + a)} L/z \\w\\ lt n + > 0. 

1 + cr 



|is(*)| < /" drllAw'CrjCfcwHIInllD^MCt) - D m i[vP]{T)\\ a 
Jo 

<£^1 + e\\i[ U i](t)\\l n 
1 + cr 



+ C(e,T)\og(l + a) f dr [\\w(t)\\\^ + ||i[u*](r) 
Jo 

<^ + eC(t) + C(e,T)log(l + a) [ £(r)dr. 
Now we consider 14(f). First we write it in the form 

h(t) = i2(t) + il(t), 

where 



(48) 



12(f) ~ / dr (D k u*3(T)D lW (T),D m i[w}(t) - D m i[w](r)) n ; 

JO 

j|(t)~ / ^(VMA^M.^HW-^HW), 
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(the star * in these formulas have the same meaning as above). Using ([47|) 
and (1421) we have 



\I2(t)\ <C T / dr\\D lW {r) (D m i[w](t) - D m i[w](r)) \\ a 
Jo 



<C T [\og{l + a)} 1 ' 2 /"*dr||«;(r)||i,n||*[ti;](t)-iH(r)|| 2i n + ^ 
Jo 

<?P- + ^(t) + C(e,T)]og(l + a) [ £(r)dr. 
1 + ° Jo 



1 + a 

(49) 



Using integration by parts we rewrite term l\{t) in the form 



Il(t) ~ / dr(u 3 (r),D k [D lW *(r) (D m i[w](t) - D m i[w}(r))}) n 
Jo 

and thus, since i? 1 / 2 (£7) C L^Tl), we have 

V b S)\< f dr\\u\r)\\ l/2 ^ 
Jo 

x \\D k [D iW *{t) (D m i[w](t) - D m i[w](T))} \\l 4/3 (U)- 

One can see that 

\\D k (fg)h 4/3{ n) < C||/IM|g||i,n, f,g £ H 1 (Jl). 

Therefore 

I4(*)l <Ct f dr\\u\T)\\ l/2 ^\\i[w](t) -iH(r)|| 2 , n 

JO 



The trace theorem implies 



KM(*)ll2,Q+ / dr\\u\T)\\{ /2M 
Jo 

' dr\\u^T)\\ 1/2 ^) + f dr\\i[w}(r)\\l n 
o J Jo 



and thus 



dr\W{r)\\i /2>n <C / dsE(i[v](s),i[v](s))) 



\I b 4 (t)\ <[e + C T , £ t] + C T , £ /" dre(r), * £ [0, T] 

JO 



(50) 



The estimates in (05J), gBJ, gSJ), (@SJ), (JgDJI and (jMD allow us to choose e 
and > such that 

£(t) < 7^- + C 2 log(l + <r) / £(r)dr for all t E [0, T#] 
l + o- Jo 

with arbitrary cr > 1. Now as in [35l [36] (see also jT2], Appendix A]) applying 
Gronwall's lemma we can conclude that £(t) = for all t E [0, TU] for some 



19 



Remark 3.5 In the case a = we can prove the existence of weak solutions 
which satisfies the energy inequality, using the same type of argument. The 
uniqueness of these solutions is still an open question. Sedenko's method 
does not work here because the nonlinearity is strongly supercritical when 
a = 0. 

Step 5. Continuity with respect to t and the energy equality. First we 
note that the vector (v(t);u(t);gul(t);guf(t);wt(t)) is weakly continuous 
in % for any weak solution (v(t);u(t)). Indeed, it follows from (f2Tj) that 
(v(t);u(t)) satisfies the relation 

(v(t),i>)o = KV(0))o + / [-vE(v,i/j) + (G f (T),i>)o]dT 

Jo 

for almost all t G [0,T] and for all ip G V = {v G V : v\n = 0} C V C V, 
where V is given by (|20p . This implies that v{t) is weakly continuous in 
Vq. Since I C l^', we can apply Lions lemma (see |31} Lemma 8.1]) and 
conclude that v(t) is weakly continuous in X. The same lemma gives us weak 
continuity of u(t) in W. Now using (|2ip again with tp G V we conclude that 

t i-> (M a w t (t), 5)n + o(u t (t),f3)n is continuous 

for every (3 = (J3; 5) G x J3j(n) x flg(n). This imply that 

1 1 — ^ (gul(t); guf(t);wt(t)) is weakly continuous 

in y = L 2 (0) x L 2 (S7) x flJ(O) for every e > 0. 

In the proof the energy equality we follow the scheme of [23]. To this 
end we need to introduce finite difference operator D^, depending on a small 
parameter h. 

Let g be a bounded function on [0, T] with values in some Hilbert space. 
We extend g(t) for all t G M by defining g(t) = g(0) for t < and g(t) = g(T) 
for t > T. With this extension we denote 

9h(f) = 9(t + h)- g(t), g^{t) = g(t) - g{t - h), 
D h g(t) = ^(g+(t)+g-(t)). 

Properties of the operator are collected in Proposition 4.3 |23| . 

Using weak continuity of weak solutions, we can extend the variational 
relation in (|17p on the class of test functions from Ct (instead of C%) by an 
appropriate limit transition. More precisely, one can show that any weak 
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solution (v;u) (with u = (u ,u ,w) = (u;w)) satisfies the relation 



L 



T i-T r-T 

(v,4> t ) dt + v\ E(v,(f>)odt- [(M a w t ,d t )n + g(ut,b t )a]dt 



10 J Jo 

+ 7/ (M a w t ,d) Q dt+ (Aw,Ad) Q dt+ a(u,b) Q dt + q(u,b) n dt 
Jo Jo Jo Jo 

= («o,0(O))o + (M«ti^,d(0))n + e(«i,6(0)) n 

- [(v(T),<f>(T)) + (M Q w t (T),d(T)) n + g(u t (T),b(T)h] 

+ / (G f (t),<f>) dt+ [ (G sh (t),b) n dt (51) 
Jo Jo 

for every <fi € £t with </>L = b = (b 1 ; b 2 ; d) = (b; d). 
Now we use 

as a test function in (|5ip . For the shell component we have test function 
& = <f>\n = -0/1^ _ the same one that used in |23j for the full Karman model. 
Thus, all the arguments for the shell component in our model are the same 
as in [23J, and we need to treat the fluid component only. Using Proposition 
4.3 [23], one can conclude, that 

t+h 

v(r)dT 

h 

l[v(T)-v(0)] + u £ dtE(v(t),v(t)) 

when h — > 0. This makes it possible to prove the energy equality in (|25|) . 

Continuity of weak solutions with respect to t stated in (|23|) and (121 
can be obtained in the standard way from the energy equality and weak 
continuity (see [HU Ch. 3] and also [23]). 

Step 6. Continuity with respect to initial data. First we prove the conti- 
nuity with respect to weak topology. 

Let g > and {Uq} be the sequence of initial data such that Uq — 1 Uo 
weakly in H (defined by (j22j) ) as n — » 00. We need to prove that 

U n (t) U(t) weakly in U for every t > 0, 

where U n {t) = {v n (t)-u n {t);v%(t)) and U{t) = {v{t);u(t);u t (t)) are weak 
solutions to the problem in question with the corresponding initial data. 
Using energy relation (|25l) and Proposition 13.41 we conclude that 

sup\\U n (t)\\ 2 H <C T , n=l,2,... 

[0,T] 



J dt (v(t),D h v(t))o + uE(v(t),^J 
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This enables us to extract a subsequence such that 



U n {t)^U*(t) * -weakly in £00(0, T; ft); 
v n {t) u*(t) weakly in L 2 (0,T;V); 

weakly in L 2 (0, T; ff Q ); 



(52) 
(53) 
(54) 



for some U*(t) = (v*(t);u*(t);u^(t)). Performing limit transition in (1171) 
and using weak convergence of initial data, we obtain that (v*;u*) is a weak 
solution and thus U*(t) = U(t) by the uniqueness statement. 

To proceed with the proof, we need to establish additional estimates for 
time derivatives. 

Taking ip G such that ip(Q) = and = as a test function in 
(fTT|) . we obtain 



which allows to define as a functional on Vq = {(f) G V : 0|n = 0} and 
get the estimate 



Using N b with b = {b 1 \b 2 -d) G C<}(0,T;iZ£(n) x fl-}(fi) x xflg(«)) as a 
test function, we obtain the estimate 

||(1 - aA)u> tt || L2 ( ,T;J-_ 2 (Q)) + Q\\Utt\\L 2 {0,T;H-l(n)) < C(T), 

where J r _ s (0) with s G [1,2] is adjoint to Hg(Q) with respect to duality 
generated by H a (Q). Thus we have the following additional convergence 
properties: 



We note that H a (Sl) C Hg(Sl) C L 2 (J2) C J r - 2 (^) for every <r G [0,1]. 
Therefore from the relations above and also from (|52p -( [54|) using the Aubin- 
Dubinsky theorem (see [37]) we conclude that 



v n -»■ u strongly in C(0, T; £T e ), 

u n -> u strongly in C(0, T; H^ e {n) x H^- £ (Q) x £$- 6 (fi)), 
< -»• w t strongly in C(0, T; F^ e (0)), 
qu? -»• gut strongly in (7(0, T; [#- £ (Sl)] 2 ). 

In particular, this implies that U n (t) — > U*(t) weakly in 7~L when Uq — > Uq. 

To prove strong continuity with respect to initial data, we use an idea 
borrowed from |23| . Due to weak continuity already established we need 





^™IIl 2 (0,T;(Vo)') ^ 



v n ^ Vf weakly in L 2 (0,T; [iJ-l(0)]3) 5 

^ ->> Qu tt weakly in L 2 (0, T; [H" 1 ^)] 2 ) 
u$ w tt weakly in L 2 (0, T; JF- 2 (H)) , 
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only to show the convergence of ||?7 n (t)||-^ to for each t > under 

the condition \\Uq — Uo\\n — > as n — > oo. 

Let Qo(u) be given by the first line of (|27|l with 

£11 = £?1 = U X1 , £22 = ^22 = u x 2 ' £ 12 = e \2 = u l 2 + 



Due to the argument given in the proof of Proposition (j3.4j) y Qo(u) is an 
equivalent norm on JT^(n) x JET^(n). Thus, to establish convergence of U n (t) 
to C/ (i) in "H it is sufficient to show that 

£ (U n (t)) -> 5 (^(*)) as n ^ oo for every f > 0, 

where 

£b(«,«,«t) = - [\\vf + \\M^ 2 w t \\l + ell«t|ln + II &M\h + Qo(u)] • (55) 



To prove this we use the energy relation in f)25[) . First we note that the 
potential energy term Q(u) in the energy £ given in ()26p has the form 

Q(u) = Q (u) + Comp(u), (56) 

where Comp (u) is a functional which is continuous with respect to weak 
convergence in %. Since £(Uq) — > £(Uq), it follows from energy equality 
(l25j) that 

lim [£ Q (U n (t)) + <b t (U n )) = £ (U(t)) + v f E(v,v)dr + $ t (U), 

where ^ t 

$ t (U) = v E(v,v)dr + j (M a w T ,w T )dT 



Using lower semicontinuity of the functional $>t(U) with respect to weak 
convergence U n (t) — ^ U(t) in 1,2(0, T;7i), we obtain that 

liminf£ (^ n (i)) < £o(U(t)), 

n— >oo 

which together with lower semicontinuity of the (quadratic) energy £q gives 
us the desired result. This completes the proof of Theorem 13.31 



4 Stationary solutions 

In this section following ideas presented in [13] and |41| we describe prop- 
erties of stationary solutions in the case when the forces Gf and G s h are 
autonomous. These solutions are the same in both cases g > and g = 0. 

It follows from Definition 13.11 that a stationary solution (v;u) G V x W 
satisfies the relation 

uE(v, if)) + (Aw, A5)n + a(u, fj) 

+ q(u,P)-(Gf,fJ>) o -(G Bh ,P) n = (57) 
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for all V G V with = p = (/3 1 ; /3 2 ; 5) and /3 = (/3 1 ; /3 2 ). The space V is 



given by (|20|) . i.e., 

V = G V | V|n = /3 = /3 2 ; 5) G W } , 

where 

w = x flj(n) x Ho(n). (58) 

Moreover, by the comparability condition we have that v\qo = 0. 
We start with the following description of stationary solutions. 

Proposition 4.1 A couple (v;u) € V x W is a stationary solution if and 
only if 

• v lies in Vq = G V : ipldO = 0} and satisfies the relation 

uE(v, %l>) - (Gf, i/>) o = 0, V V G V i (59) 

• u = (u l ;u 2 ;w) £W = Hq(Q) x Hq(Q) x Hq(Q) satisfies the equality 

(Aw, A5) n + q(u, (3) - (G sh + N^G f , 0) a = (60) 

for all (3 = (/3 1 ; (3 2 ; (5) G W, w/iere W is given by ([58)) and 

g(«,/3) =(^11 + ^22,^ (61) 
+ (Nnw xi + N 12 w X2 ,S xl )n + (N 12 w Xl + N 22 w X2 ,5 X2 )n 

+ (iVn,/4)n + (JVu,^ +/5^)n + (iW 2 2 )n 



Proof. Taking /3 = in (|57|) yields (|59j) . Since v\ QO = 0, one can see that 

E(v, Nof3) = for every f3 € W. Therefore taking ip = Nq(3 in (|57l) gives us 
(|60p . Thus any stationary solution satisfies (|59p and (|60p . Similarly, we can 
derive (EH) from (ED and d6QD. □ 



The following assertion shows that for the forces Gf and G s h of a special 
structure we can guarantee the existence of stationary solutions. 

Proposition 4.2 We assume that 

G f = 0, Gl h = G 2 sh = and G 3 sh = g G H~ l (O) . (62) 

TTien any stationary solution has the form (0; u), where u = (u ; u 2 ; w) G 
satisfies 

(A™,A£) n + ?(u,/3)-(<7,£) n = 0, V/3= (/^V) G W. (63) 

Moreover, there exists at least one solution u = (u ,,u 2 ;w) to (j63|) in the 
space W . If in addition we assume that k\ = k 2 = 0, then the set of all 
stationary solutions to (|63p from W is bounded in the space W . 
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Proof. If Gf = 0, then it follows from (|59|) that v = 0. Thus equation 
(1601) turns into ([63]) . 

To prove the existence of elements u = {u l ;u 2 ;w) £ satisfying (|63p . 
we note (see, e. g. [H]) that a solution can be obtained as a minimum point 
of the functional 

U(u) = ^[\\Aw\\ 2 + Q(u)]-(w,g)n on 

where Q is given by (|27p . It is clear that n(u) is bounded from belowH 
Thus we can construct appropriate Galerkin approximations {u n } C W for 
the global minimum and note that Tl(u n ) < 11(0). This facts together with 
Proposition 13.41 provide us an a priori estimate which allows to prove the 
existence of a solution by the same method as in [41] . 

To prove the boundedness of stationary solutions we consider the func- 
tional q(u,(3) given by (itITj) with j3 = (/3 1 ; /3 2 ; 5), where with 5 = \w and 
f3 l = u l . In this case we obtain that 

q(u,p) =- hkiNu + k 2 N 22 ,w) n + (N ll ,u 1 xl + hw + ^[w Xl ] 2 ) n 

+ (N 12 ,u X2 + u 2 Xl + w x1 w X2 )q + (N 22 ,u 2 X2 + k 2 w + ^[w X2 ] 2 )q. 
Using the expressions for the deformation tensor {£ij} we obtain 

q(u,f3) =- -(kiN n + k 2 N 2 2,w)n 

+ (N n ,e n )n + (Ni 2 ,ei 2 )n + (N 22 , e 22 ) n . 
The Hooke law (|4|) yields 

q(u, /?) = - -(fcliVn + k 2 N 22 , W ) Q 

+ (1 -/x) j n H 1 + ^2 + 2 ^ £ H £ 22 + ^(1 - M)£l 2 ) 

= - -(hN n + k 2 N 22 , W ) Q + Q(u), 

where Q(u) is given by (|27p. Therefore under the conditions in ([62p from 
([63]) we have that 

-||Au>||n + Q(«) " ^(^1^11 + A; 2 iV22,u»)n = \{g,w) n 

Thus, the set of stationary solutions is bounded provided ki = (or even 
small enough). □ 

To the best of our knowledge the existence of stationary solutions in the 
case of general external loads is still an open question. 



2 This is exactly the point where we use the structure of the external forces assumed in 
(ED 



25 



5 Existence of a global attractor 



In this section we prove the existence of a compact global attractor under 
the condition that the external forces satisfy (|62|) . 

First we note that by Theorem 13.31 problem jH)-© generates an evolu- 
tion semigroup St in the space 7i, which has the form 

• H = X x W x Y in the case g > 0, 

• H = X x W x H a in the case g = 0, 

where X and Y are defined in ([lip and (|12j) and W is given by (|58|) . The 
evolution operator St is defined as follows 

• Case g > 0: St(«o; i*o; u{) = U(t) = (v(t);u(t);ut(t)), where the 
couple (v(t);u(t)) solves (HHP- 

• Case £ = 0: St (vq ; uq ; w i ) = U(t) = (v(t);u(t);wt(t)), where v(t) and 

= (u 1 ^);^); «;(*)) solves ©-© with g = 0. 

Our main result in this section is the following theorem. 

Theorem 5.1 Assume that 7 > and the external forces satisfy (|62|) , Let 
the set of the stationary points in ~H of the problem ([I])-® is bounded. Then 
the evolution semigroup St generated by this problem possesses a compact 
global attractor. 

We recall (see, e.g., [Hll0[l38]) that global attractor of the dynamical system 
(St, rl) is defined as a bounded closed set 21 C % which is invariant (St 21 = 21 
for all t > 0) and uniformly attracts all other bounded sets: 

lim sup{dist^ (Sty, 21) : y € B\ = for any bounded set B in T~L. 

Proof. It follows from energy inequality in (j25j) that the set 

Wr = {U : $([/) = £([/) - (<7,w)n < i?} 

is forward invariant with respect to St for each R > 0. Here ?7 = (v;u;ut) 
with n = (li 1 ;^ 2 ;^;) in the case g > and f7 = (v;u;wt) in the case £> = 0. 
As in the proof of Proposition 14.21 using (162ft one can see that &(U n ) — > +00 
if and only if ||[/ n ||% — > +00. Therefore the set Wr is bounded and any 
bounded set belongs to Wr for some R. Moreover, it follows from energy 
inequality (f25j) that the continuous functional $(U) on H possesses the 
properties (i) ®(S t U) < $(U) for all t > and U G H; (ii) the equality 
<&([/) = $>(SfU) holds for all f > only if U is a stationary point of St. This 
means that 3>(J7) is a strict Lyapunov function and (7i,St) is a gradient 
dynamical system. Therefore due to Corollary 2.29JTT] (see also Theorem 4.6 
in [M]) we need only to prove asymptotic smoothness. We recall that (see, 
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e.g., [22J or [TT]) that a dynamical system (X, S t ) is said to be asymptotically 
smooth if for any closed bounded set B C X that is positively invariant 
(StB C B) one can find a compact set fC = fC(B) which uniformly attracts 
B, i. e. sup{distx (Sty, /C) : y € B} — > as t — > oo. 

To prove asymptotic smoothness we use Ball's method of energy relations 
(see [5] and also |32j). For a convenience we recall the abstract theorem (in 
a slightly relaxed form) from |32j which represents the main idea of the 
method. 

Theorem 5.2 (|32j) Let St be a semigroup of strongly continuous opera- 
tors in some Hilbert space rl. Assume that operators St are also weakly 
continuous in rl and there exist a number uj > and junctionals A, L and 
K on rl such that the equality 

A(t) + f L(T)e- 2w ^dT = A{s)e~ 2 ^ t ~ s) + f K^e-^-^dr, (64) 

J s J s 

holds on the trajectories of the system (H,St). Here we use the notation 
G(t) = G(StU), where G is one of the symbols A, L and K . 
Let the functionals possess the properties: 

(i) A : TL — > M + is a continuous bounded functional and if {Uj}j is 
bounded sequence in % and tj — > +oo is such that (a) St 3 Uj — U 
weakly in H, and (b) limsup n _ i>00 A(St j Uj) < A(U), then S^Uj — > U 
strongly in T-L. 

(ii) K : % — » R is 'asymptotically weakly continuous' in the sense that 
if {Uj}j is bounded in %, and St^Uj — 1 U weakly in % as tj — > +oo, 
then K(S T U)£ L[ oc {R + ) and 



lim I e- 2uj( - t -^K(S s+t U j )ds= f e~ 2uj{t -^K(S s U)ds, Vi > 0. 
Jo Jo 

(65) 



(iii) L is 'asymptotically weakly lower semicontinuous' in the sense that 
if {Uj}j is bounded in %, tj — > +oo, SuUj — 1 U weakly in H, then 
L(S T U) G L' 1 oc (E + ) and 

liminf f e-M*-') • L(S s+t .Uj)ds > [ e - 2u(t ~ s) L{S s U)ds, Vt > 0. 
Jo Jo 

Then St is asymptotically smooth. 

Now we check validity of the hypotheses of Theorem 15.21 in our case. 

Step 1 (continuity): The continuity of the evolutional operator in both 
(strong and weak) senses follows from Theorem 13.31 

Step 2 (energy equality): All the calculations below can be justified by con- 
sidering approximate solutions. 
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Let = ^(v(t),u(t),u t (t)), where 

i£(v,u,ut) = (M a w t ,w)si + g(u t ,u) n + («,iVo[u]) , 

where f7 = (v; u; ut) = (v; u; w; u t ;wt) is a weak solution and Nq is given by 
(fT6j) . One can see that 

jj®(t) =(M a w t t,w)n + g(u a ,u)n + (ut,iVoM)o 

+ ||MyV||E + *lln + (^^ok])o 
=||Mi /2 ™t|ln + e \\u t \\l + {v,N [u t ])o - uE(v,N [u]) 

- \\Aw\\ 2 - Q(u) -j^—f [< + < + (1 + ^' 

-j(M a w t ,w) n + <S>o(t) + (g,w) Q , 
where &o(t) is a linear combination of the terms of the form 

(Diu l DjW,D m w)n and (kiwDiw, D m w)n- 
Let A(t) = £{t) + r]^(t), where r] > will be chosen later. Since 

j t E{t) = -uE(v,v) - iWM^wtWl + (g,w t ) n 



one can see 



j t A(t) + (7 - v)\\M^ 2 w t \\ 2 n + vE(v,v) - riQWutWh 



+ 1] 



||Aw;|| 2 + Q(«) + — !— / [«4 + w 4 X2 + (1 + ^)«] ds ' 
1 _ I 1 Jn 



where 



=77[(v,iVoH)o-i/E(u,iVo[n])-7(M a ^,K;)n] 
+ »7*o(*) + ^(3, w)n + (5, ^)n. 

Consequently 

j t K + 2wA + (7 - 77 - ^llMyVll 2 } + (1/ - u) - (r? + u) g \\u t \ 

+ { V -u) [||Au;|| 2 + Q(n)] 

+ -r— [ Ki + < + C 1 + /*)«] ^ = + ^(t). 
1 - V Jn 

Thus we obtain (|64p with L and if given by 

=( 7 _ rj - u:)\\M l J 2 Wt \\l + oo)E(v,v) 
- ( v + oj)g\\u t \\ 2 n + ( V -u) [\\Aw\\ 2 + Q(u)] 

Tj J 

+ 1 / [*>xi + w t 2 + (1 + V) w2 xi w l 2 ] dx> , 

l ~ V Jn 
K(t) =^ 1 {t)+urj^(t). 
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We choose 7] > oj > such that 7 — 7] — uj > and v — to > 0. 

Step 3 (properties of the junctionals): Now we prove that the functionals A, 
L and K satisfy requirements (i)-(iii) in Theorem 15.21 First we rewrite the 
energy £ in the form 

£(v, u, u t ) = 6q(v, u, u t ) + Comp (u), 

where the quadratic energy functional £q is given by ([55)) and Comp (u) 
is a functional which is continuous with respect to weak convergence. By 
Proposition 13.41 the functional £q(v, u, ut) provides an equivalent norm on % 
and therefore £{t) satisfies (i) by the properties of weak convergence. Now 
we show that \P is weakly continuous. We start with the third term. If 
Uj — 1 u weakly in W, then due to Proposition 12.21 NqUj — Nqu weakly in 
[H 3 / 2 (0)} 3 f]X and strongly in X. Thus, (v j: N Q v?)o -»• (v,N u) O - The 
remaining terms are obviously weak continuous. Thus A(t) satisfies (i). 

Now we consider K. As above ^ and all the terms in $1 (except of 
$0) are obviously weak continuous. The same is true for $0 due to the fact 
that (/i;/2) /1 • /2 is a (strongly) continuous from i? 1//2 ($7) x // 1 / 2 (0) 
into L<i{£l). Applying this property to the terms (Diu l DjW, D m w)n and 
(kiwDjw, D m w)Q, we find that the functional K is weakly continuous and 
thus by the Lebesgue dominated convergence theorem the convergence of 
integrals in (I65p holds. 

As for the functional L, all its terms are obviously weakly lower semi- 
continuous (because of the convexity norm properties, relation (I56h and 
Proposition I3.4p . except of E(y, v) and — ^H - ^ ||^. Let {Uj} C H is bounded, 
tj — > 00, and StjUj — weakly in %. Denote by Uj(s) a solution to 
the system under consideration with the initial data St Uj. From the en- 
ergy balance equality in ([25]) we have that the sequence of the velocity field 
components Vj of Uj(s) is bounded in ^(0, t;V). Therefore due to weak 
continuity property of the form E(v,v), we obtain 

f e- u{t ~^E{v{s), v{s))ds < liminf I e^-^ E(vj{s),Vj(s))ds. 
Jo Jo 

Due to the standard trace theorem we can suppose that dtUj = u||n) 

1/2 

converges to dtu weakly in Hj (f2)) (hence strongly in L2(fi)) for almost 
all t > 0. Therefore the Lebesgue theorem yields that dtUj — > dtu strongly 
in £2(0, T; £2(0)), so t ne property (hi) holds. 

Thus we can apply Theorem 15.21 to prove asymptotic smoothness and 
complete the proof of Theorem 15.11 □ 

Remark 5.3 We do not know whether we can relax essentially the struc- 
tural force hypothesis ([62]) in Theorem 15. 11 We use ([62]) to prove the bound- 
edness of the sublevel set Wr only. The question on the validity of this weak 
form of the coercivity of the energy functional is still open in the case of 
general external loads. 
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